Stabilization of Membrane Pores by Packing 
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We present a model for pore stabilization in membranes without surface tension. Whereas an 
isolated pore is always unstable (since it either shrinks tending to re-seal or grows without bound 
til to membrane disintegration), it is shown that excluded volume interactions in a system of many 
pores can stabilize individual pores of a given size in a certain range of model parameters. For such a 
multipore membrane system, the distribution of pore size and associated pore lifetime are calculated 
within the mean field approximation. We predict that, above the temperature Tm when the effective 
line tension becomes negative, the membrane exhibits a dynamic sieve-like porous structure. 



Pores can form and grow in membranes in response 
to thermal fluctuations and external influences. Pore 
growth enhances the transport of biomolecules across 
the membranes and its biological relevance can bring 
new prospective biotechnological applications (see, e.g., 
Rcf. 0, I2I). Holes appear in the membrane via a 
thermally activated poration process, and their subse- 
quent growth is controlled by the effective line ten- 
sion (assuming a negligible small surface tension). A 
schematic depiction of a membrane with holes is dis- 
played in Fig^ Preliminary computer simulations of 
coarse-grained models for lipid membranes have indi- 
cated that close to their disintegration point, membranes 
sometimes exhibit a perforated state where they are pep- 
pered of many pores Q- In the present paper, we inves- 
tigate a mechanism that can stabilize such a state. 

We first consider a single pore formation in a mem- 
brane with zero surface tension. For simplicity, the 
pore edge is regarded as a closed self-avoiding path of 
n steps of constant length I. In the simplest approxi- 
mation, the free energy of the system (tensionless mem- 
brane with a pore of size n) at nonzero temperature 
can be written as a summation of two contributions: a 
purely energetic part as suggested by Litster [J|, plus 
an entropic part as modeled by Shillcock and Boal Q, 
/(n) = fo + XqIti — kBTln[Lj(n)], where /o is a n- 
indepcndent energy, Aq the bare line tension of the pore 
edge, ksT = l/f3 the thermal energy, and uj{n) the num- 
ber of possible conformations of the pore contour of size 
n. For self-avoiding walks, Lu{n) has a general form [3,13 
given by uj{n) = ujq z" 71"^^, where ujq is a constant, z the 
connectivity constant of the medium, and the exponent 
a = 1/2 for self-avoiding random walks in two dimen- 
sions and a < 2 for any kind of pores. Let no be the 
minimal size of the pore, the above free energy can be 
rewritten as f{n) = for < n < no, and for n > uq, 



f{n) = Fo + Xo (l-^] l[n~no] 



(2-a)kBTln — 



(1) 



where we have defined Fq as the free energy required 
to create or initiate a minimal-sized pore, Ai the en- 
tropically modified line tension of the pore edge, and 
Tm = Ao V^B lii(z) the disintegration temperature. In 
what follows, we will neglect the logarithmic term, (2 — 
a)ln(n/no), that only slightly rcnormalizes the results. 
Simple inspection of Eq.lQJ indicates that the free en- 
ergy of the system monotonically increases as the pore 
grows larger at low T < Ti„ when Ai > 0, keeping the 
membrane stable with an unstable pore that reseals. In 
this case, a pore of any size has a finite lifetime and will 
ultimately shrinks to disappearance. In contrast, when 
the effective line tension is negative, Ai < 0, at high 
T > Tm, all newly initiated pores grow without bound 
(i.e.; with a diverging lifetime) leading hence to destabi- 
lization and disintegration of the membrane. There are 
numbers of works that have investigated mechanisms for 
stabilizing membranes with a single pore. These include, 
for instance, membrane bending fluctuations, renormal- 
ization of linear and surface tension coefficients [g , area 
exchange in tense membranes |9j, osmotic stress |lOJ . 
hydrodynamics [Hi Il2l | , orientational ordering [l3J , and 
others (sec Ref. J14j for more details and references). In 
this paper we focus on a new stabilizing effect, which is 
created by the presence of multiple pores in the mem- 
brane. 

Indeed, as the probability of initiating several pores 
on a membrane increases as the temperature gets higher, 
we are now confronted to the situation of a membrane 
with an ensemble of pores (as in Ref. [ia|)- If the bare 
line tension is negative, the system tries to create as 
much pore rim as possible |l6l |. In that case, a mem- 
brane state with many small pores is more favorable than 
one with only one very large pore. In some sense, such 
a state is similar to a droplet microemulsion structure 
in amphiphilic systems, where the fluid is macroscop- 
ically homogeneous, but filled with internal interfaces 
on the microscopic scale [l7l |. To quantify this expec- 
tation, we proceed as follows: We first consider the free 
energy F{N{n)} of a membrane where the number N{n) 
of pores with contour length n is fixed (constrained). In 
the second step, we will relieve the constraint and min- 
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FIG. 1: Schematic depiction of a porous membrane, illus- 
trating the different areas introduced in the text: The total 
surface Aup covered by all lipids (gray filled circles), the total 
area of pores Aporos, and the total area Am ~ ^iip + j4poros- 
The circle around the pore represents the area that is effec- 
tively excluded by the pore, and the total excluded area for all 
pores is Aox- Due to the fractal shape of pores, the excluded 
area of pores is larger than their actual area. 



imize F{N{n)} with respect to N{n). The free energy 
F{N(n)} has energetic line tension contributions and en- 
tropic shape contributions as in Eq. (^. Moreover, the 
pores have the translational entropy of a two dimensional 
gas. However, they may not overlap, since a configura- 
tion with two "overlapping" pores would have to be re- 
placed by a new configuration with just one, larger pore, 
and a different N{n). For fixed N{n), this restricts the 
translational degree of freedom of pores, as if they had 
excluded volume interactions. Note that these "interac- 
tions" are purely entropic. One could also introduce real 
repulsive interactions between pores, originating, e.g., 
from direct electrostatic or van der Waals forces, entropic 
or Helfrich interactions arising from pore shape fluctua- 
tions, or from membrane undulations |l8lll9l|. However, 
this is not necessary for our argument. 

To proceed we adopt a simple van der Waals ap- 
proach and approximate the free energy F{N{n)} by 
F{N{n)} = En^HfH - keT In^, where f{n) is 
the free energy of a single pore as given in Eq.J^I, and 
the total translational entropy fl of the pores is given by 



ii [N{n)]\ 
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Here ao is an area constant defined below, A^ is the to- 



tal membrane area, A,^ 



Aiip + Aporcs, with the area 
and Apx is the 
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of lipids A\ip and the pore area A 
area that is effectively inaccessible to a test pore due 
to the presence of the other pores (the excluded area). 
Due to the fractal nature of the pores, ^ox is larger than 



Spores as indicated in Fig. ^ Notice that A\ip is con- 
stant since it is proportional to the number of lipids. 
In writing the rightmost expression in Eq.|(21), we have 
used the approximation, A^! « N^ e~^ , and rewritten 
the accessible area for pores as A,„ — A^^ = Aup — Aq, 

such that Aq = Aox — Spores = y^N{n)a{n), where 

n 

a{n) is the difference between the excluded and the ac- 
tual areas of a pore of size n. Since the contours of the 
pores have self-avoiding walk statistics, a{n) scales like 
a{n) = ao(n/no)^'^ with the Flory exponent i/ « 3/4. 
This defines oq. Now, inserting Eq.Q into the free en- 
ergy expression, and minimizing with respect to N{n), 
yields the normalized equilibrium distribution Pcq(^) of 
pore sizes. 



Peq(n) = iV(n) /^iV(n) 



exp{-/3G(n)} 

Q 



(3) 



Here we have defined the effective free energy, G{n) = 
f{n) -\- 'k-Q'TQin / nof'^ . The partition function Q can be 
regarded as a packing density of pores, and is given self- 
consistently as. 



ao 
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-^^^^Ux, (4) 



Defining x ~ n/no and using i/ = 3/4, the effective free 
energy G{x) of a pore in a membrane containing an en- 
semble of pores (in gas phase) reads as G{x) = for 
< X < 1, and, 

G{x) =Fo + QkeT + nol f Xoff{z) dz ; x>l, (5) 

where the effective line tension Aeff is defined as, 



Aeff(a;) =: Ao 1 



T 



/SQkeT 



.1/2 



(6) 
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As a result of the renormalization of Eq.(^ by the ex- 
cluded area constraints due to other pores, the effective 
line tension becomes a function of the pore size. This 
is the origin of the packing stabilization mechanism that 
leads to sieve-like structure of membranes. Indeed, as 
already mentioned above, a single pore destabilizes the 
membrane at high temperature T > T^ as the effective 
line tension becomes negative. In a multiple pore sys- 
tem, however, the membrane remains stable even beyond 
r,„. For temperatures such that Tm < Tc < T, G{x) in 
Eq. |SJ) admits a minimum at x = xi such that 



Aoff(a;i 
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T-Tn, 

T, - T„ 



(7) 



where Qc = Q(T = Tc) and Tc, the critical temperature 



at which xi coincides with x 
2Qo(l - Tn,/Tc), i. e., 

2Qo(i-y) 



1, is given by 3Qc 



-yPmFo 



(8) 
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FIG. 2: Reduced equilibriuin distribution, Peq(a;)/Pcq(l) ~ 
g-/3[G(r.)-G(i)l^ of pore size x for /3Fo = 2, Qo = 4 and two 
temperature regimes. G{x) is given in Eq.0 with Q — 0.215 
for T = T™ and Q = 8/(9^/3) = 0.513 for T > T^ (i.e., for 
T = 3r™^/2 and Tc = 1.082T„). The dashed vertical line 
indicates the maximum xi = 3 for T > T^. 



with Qo = '^o ln(z) and y = Ti„/Tc. Interestingly, Eq.Q 
states that the sizes corresponding to a vanishing line 
tension is that of stable pores. 

Two conclusions can be drawn from these results. 
First, the excluded area constrains between pores can 
stabilize a membrane even in parameter regions where 
the effective line tension of pores is negative. Second, the 
distribution of pore size in the case of negative effective 
membrane line tension may have a maximum at nonzero 
contour length as illustrated in Fig. [21 Three regimes can 
be distinguished: (i) At low T < r,„, the line tension is 
positive, and the distribution Pcq{n) of pore size drops 
monotonically as a function of n. (ii) At intermediate 
temperatures, T„i < T < Tc, the line tension is nega- 
tive, Pcq{n) still drops monotonically, but the pores are 
now stabilized due to the presence of the others, (iii) At 
high T > Tc, a maximum emerges in Peq(n), i. e., pores 
have a most probable size. Figure |21 shows the reduced 
distribution of pore size, Pcq{x) / Pcq{^) , for two temper- 
atures below and above Tc, with a maximum at x = 3 (or 
n = 3no) for T > Tc- This regime sensitively depends on 
the bare pore free energy, Fq, required for creation of a 
minimal-sized pore. This is illustrated in Fig. |31 where Tc 
decreases towards T^ (i.e., Tm/Tc increases towards 1) 
when increasing either Fq or the effective smallest pore 
size Qo- 

From a practical point of view, the porous membrane 
state is interesting because it should have peculiar perme- 
ability characteristics. Apart from the number and size 
of pores, another quantity determining the permeability 
is the lifetime of pores, i. e., how long a pore stays open 
once created. We now discuss briefly how the lifetime of 
pores is affected by the presence of the other pores. 

Six processes contribute to the dynamical evolution of 
the pores: Pore opening and closing, pore growth and 
shrinking, pore coalescence and splitting. The pore open- 
ing and closing is mainly controlled by the potential bar- 
rier that must be overcome to create a pore: the am- 
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FIG. 3: Ratio Tm/Tc in Eq.Q as a function of the reduced 
energy of pore formation. 



phiphiles must change their orientation, and the free en- 
ergy in the intermediate state is different from the corre- 
sponding energy when a pore already exists. The char- 
acteristic time scale of this process does not depend on 
the pore packing and shall not be considered here. Am- 
phiphile rearrangements also take place when pores coa- 
lesce and split, hence these processes are rate-driven and 
the characteristic time scale depends on the height of a 
potential barrier. For simplicity, we shall assume that 
this height is very large, i. e., we shall neglect coales- 
cence and splitting events. Hence we arc left with the 
characteristic time r for the growing and shrinking of a 
pore. Our main interest in this quantity is to allow com- 
parison on the pore lifetime in two distinct situations: a 
single pore in a membrane versus a pore in a gas phase 
ensemble of pores in a membrane. 

Neglecting hydrodynamic effects J23, the growth and 
shrinking dynamics of the pore can fairly be described by 
a diffusion process with a diffusion constant D (assumed 
here independent of x) in the potential G{x). Then, r 
is given by the mean time that an already existing pore 
takes to first reach the minimum pore size x = 1 (or, no). 
According to the first passage time theory the lifetime r 
can be estimated from the relation l2lL 122. Hal 
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In the case of a single pore in a membrane, Eq.(|5J) reduces 
to the pore lifetime tq as. 



To 



D 



Qo 



T 



1 



; T<T„ 



(10) 



This To quadratically increases with temperature and di- 
verges with the membrane disintegration as T approaches 
Tm. As displayed in Fig. 01 in the case of an ensemble of 
pores, the pore lifetime r begins to increase considerably 
at Tjn but without diverging and, thus with no membrane 
disintegration. At T below T^,, a pore of any size will ul- 
timately shrink to disappearance on average because of 
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FIG. 4: Reduced pore shrinking time in Eq.Q as a func- 
tion of tlie reduced temperature (3o(l — Tm/T) for different 
/3-Fo- The dotted vertical line indicates T = Tm above which 
the line tension becomes negative. Dashed vertical lines in- 
dicate the crossover points T = Tc [i.e., Qo{l — Tm/Tc) — 
0.975, 0.771, 0.463 for pFo = 0.5, 1, 2, respectively] where 
Pcq{x) starts to exhibit a maximum like in Fig. |^ 



the drift towards a; — 1. In the stabihzed regime above 
Tci the pore size diffuses towards the minimuni xi > 1 
of G{x), and x = \ can be reached via an escape process 
over the energy barrier, A ~ G(l) — G{xi), such that r 
scales as T oc exp(/? A). As a result, the pore remains 
open much longer in the packing stabilized regime with- 
out membrane disintegration. 

To summarize, we have analyzed the statistics of mul- 
tiple pore formation in membranes without surface ten- 



sion within a classical van-der Waals approach. We 
have shown that the exclude area interactions between 
pores allows membrane stabilization against disintegra- 
tion even when the effective line tension becomes negative 
at high temperatures. In a certain range of parameters, 
this leads to a nano-porous membrane state where the 
membrane has a sieve structure with long-lived holes of fi- 
nite size. Owing to the diversity of biological systems (all 
the more artificial lipid bilayers) and wide range of acces- 
sible parameters, we expect such stable multipore mem- 
branes exist at physiological conditions with pore sizes 
in the range of (1 — 10) nm. This invites speculations on 
possible applications of such structures, e.g., membranes 
with selective permeability for controlled drug delivery, 
or to promote biomolecule translocation. 

Our result illustrates the rich diversity of membrane 
structures that can form by chemically or physically tun- 
ing the line tension A. As a crude approximation, we have 
A ex 7o/i, where 70 is the free energy per unit area be- 
tween coexisting regions of hydrophobic and hydrophilic 
molecules, and h the membrane thickness. Such an esti- 
mation is valid only for the so-called hydrophobic pores, 
while for hydrophilic pores, the scaling is A ex k//i, where 
K ex /i^ is the bending modulus and 1/h the membrane 
pore curvature. As A ex /i in any case, one would ex- 
pect sieve-like membrane structures to be more likely in 
thin membranes although other mechanisms do exist, for 
example, in mixed membranes or membranes with addi- 
tives, that lead to reduction of the line tension. 
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